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Abstract. We define the notion of a finite-time singularity of a vector field and then discuss 
a technique suitable for the asymptotic analysis of vector fields and their integral curves in the 
neighborhood of such a singularity. Having in mind the application of this method to cosmology, 
we also provide an analysis of the time singularities of an isotropic universe filled with a perfect 
fluid in general relativity. 



1. Introduction 

There are two approaches to characterizing spacetime singularities in a cosmological context. 
The first approach may be called geometric and consists of finding sufficient and/or necessary 
conditions for singularity formation, or absence, independently of any specific solution of the 
field equations under general conditions on the matter fields. Methods of this sort include 
those based on an analysis of geodesic congruences in spacetime and lead to the well known 
singularity theorems, cf. pQ, as well as those which are depend on an analysis of the geodesic 
equations themselves and lead to completeness theorems such as those expounded in [2], and 
the classification of singularities in [3] . 

The second approach to the singularity problem can be termed dynamical and refers to 
characterizing cosmological singularities in a geometric theory of gravity by analysing the 
dynamical field equations of the theory It uses methods from the theory of dynamical systems 
and can be global, referring to the asymptotic behaviour of the system of field equations for 
large times, or local, giving the behaviour of the field components in a small neighborhood of 
the finite-time singularity. 

In this latter spirit, we present here a local method for the characterization of the asymptotic 
properties of solutions to the field equations of a given theory of gravity in the neighborhood of 
the spacetime singularity 1 . We are interested in providing an asymptotic form for the solution 
near singularities of the gravitational field and understanding all possible dominant features 
of the field as we approach the singularity. We call this approach the method of asymptotic 
splittings. 

1 Similar approaches have been used in the past to test given systems of equations for integrability under the 
so-called Painleve test. However, our approach is more geometric in nature and is not related to integrability. 



In the following sections, we give an outline of the method of asymptotic splittings with a 
view to its eventual application to cosmological spacetimes in different theories of gravity. For 
the sake of illustration, in the last Section we analyze the asymptotic behaviour of a Friedmann- 
Robertson- Walker (FRW) universe filled with perfect fluid in Einstein's general relativity, which 
provides the simplest, nontrivial cosmological system. 

2. Finite-time singularities: Definitions 

It is advantageous to work on any differentiable manifold M. n , although for specific applications 
we restrict attention to open subsets or R n We shall use interchangeably the terms vector field 
f : M. n — ► TAi n and dynamical system defined by f on Ai n , x = f (x), with (•) = d/dt. Also, we 
will use the terms 'integral curve' x(t,xo) of the vector field f with initial condition x(0) = xo, 
and 'solution' of the associated dynamical system x(i;xo) passing through the point xo, with 
identical meanings. 

Given a vector field f on the n-dimensional manifold M. n , we define the notion of a general 
solution of the associated dynamical system as a solution that depends on n arbitrary constants 
of integration, x(t; C), C = (c±, ■ ■ ■ , c n ). These constants are uniquely determined by the initial 
conditions in the sense that to each xo we can always find a Co = (cio, • • • , c„o) such that the 
solution x(t; Co) is the unique solution passing through the point xo. Therefore, a property 
holds independently of the initial conditions if and only if it is a property of a general solution 
of the system. 

A particular solution of the dynamical system is any solution obtained from the general 
solution by assigning specific values to at least one of the arbitrary constants. The particular 
solutions containing k arbitrary constants can be viewed as describing the evolution in time of 
sets of initial conditions of dimension k strictly smaller than n. A particular solution is called 
an exact solution of the dynamical system when k = 2 . Thus, in our terminology, a particular 
solution is a more general object than any exact solution, the latter having the property that 
all arbitrary constants have been given specific values. The hierarchy: exact (no arbitrary 
constants) to particular (strictly less than maximum number of arbitrary constants) to general 
solutions, will play an important role in what follows. 

General, particular, or exact solutions of dynamical systems can develop finite-time 
singularities; that is, instances where a solution x(t;c±, • • • , c^), k < n, misbehaves at a finite 
value t* of the time t. This is made precise as follows. We say that the system x = f(x) 
(equivalently, the vector field f) has a finite-time singularity if there exists a t* € R and a 
xq G M n such that for all M 6 R there exists an S > such that 



for \t — t*\ < 5. Here x : (0,6) — > M n , xo = x(io) for some to € (0,6), and || • \\lv is any LP 
norm. We say that the vector field has a future (resp. past) singularity if t* > to (resp. t* < to). 
Note also, that to is an arbitrary point in the domain (0,6) and may be taken to mean 'now'. 
Alternatively, we may set r = t — t*, r G (— t*, 6 — t*), and consider the solution in terms of the 
new time variable r, x(r;xo), with a finite-time singularity at r = 0. We see that for a vector 
field to have a finite-time singularity there must be at least one integral curve passing through 
the point xq of M. n such that at least one of its L p norms diverges at t = t*. We write 



2 There are solutions, called singular solutions, which have the property that a certain jacobian vanishes on them, 
and so they are not obtainable from the general solution of the system like the particular solutions. Below we 
shall use this term but with a totally different meaning as our notion of a solution with a time singularity (which 
can be also called a singular solution) is completely different. 



x(t;x )|| L P > M, 



(1) 



lim ||x(t;x )||LP = oo, 



(2) 



to denote a finite-time singularity at i*. 

One of the most interesting problems in the theory of singularities of vector fields is to find 
the structure of the set of points xo in M n such that, when evolved through the dynamical 
system defined by the vector field, the integral curve of f passing through a point in that set 
satisfies property (J2J). 

Another important question, of special interest in relativistic cosmology, is to discover the 
precise relation between the finite-time singularities of vector fields that arise as reductions of 
the field equations and those that emerge in the form of geodesic incompleteness. The difficulty 
here is that the finite-time singularities of vector fields appear to be unconnected to geodesic 
incompleteness and conversely, singularities which arise through the formation of conjugate 
points do not seem to demand, or require, any dynamical description. 

3. Fixed and movable singularities 

Finite-time singularities of (general or particular) solutions of linear dynamical systems are 
located at the singularities of their coefficients and are fixed because they are known from the 
singularities of the coefficients of the system. The fixed singularities in a solution are therefore 
independent of the choice of initial conditions. In contrast, solutions of nonlinear systems can 
develop finite-time singularities that are either fixed or movable. A singularity is fixed if it is a 
singularity of x(t; C) for all C; otherwise, we say it is a movable singularity. 

Note that any fixed finite-time singularity of a particular solution (k < n) cannot be a fixed 
singularity of a general solution since at least one of the constants appearing in the general 
solution has been set to zero, and so this singularity is not one of a x(i; C) for all C, that is 
independent of the initial conditions. Hence, fixed finite-time singularities in a general solution 
cannot be understood by studying fixed singularities in particular solutions. However, movable 
singularities of a particular solution, if they are singularities (in the sense of the definition above) 
of the general solution, will always be movable ones. Therefore, movable finite-time singularities 
in particular solutions make a nonzero contribution to the singularity pattern of the vector field 
and must be taken into consideration in the general study of its singularities. 

It may also happen that a dynamical system has no movable singularities in the general 
solution but still has a singular or particular solution with a movable singularity. Hence, a fixed 
(or movable) singularity in a general solution can be a fixed (or movable) singularity of some 
particular solution but not vice versa. In general, movable singularities are more interesting, 
since the issue of choosing initial conditions plays an important role for them; consequently, we 
shall restrict our attention to them almost exclusively in what follows. 

How should we tackle the geometric problem of describing the behaviour of vector fields and 
their integral curves - solutions of the associated dynamical system - in the neighborhood of a 
finite-time movable singularity? Assume that we are given a vector field, and we know that at 
some point, t*, a system of integral curves, corresponding to a particular or a general solution, 
has a (future or past) finite-time singularity in the sense of definition (fT]l. The approach we 
take in this paper is an asymptotic one. The vector field (or its integral curves) can basically 
do two things sufficiently close to the finite-time singularity, namely, it can either show some 
dominant feature or not. In the latter case, the integral curves can 'spiral' in some way around 
the singularity ad infinitum so that Q is satisfied and the dynamics are totally controlled by 
the subdominant (lower-order in terms of weight - see below) terms, whereas in the former case 
solutions share a distinctly dominant behaviour on approach to the singularity at t* determined 
by the most nonlinear terms. 

To describe both cases invariably, we decompose the vector field into simpler, component 
vector fields and examine whether the most nonlinear one of these shows a dominant behaviour 
while the rest become subdominant in some exact sense. Using this picture, we then built a 
system of integral curves corresponding, where feasible, to the general solution, and sharing 



exactly its characteristics in a sufficiently small neighborhood of the finite-time singularity. The 
construction of these solutions is given as a formal asymptotic series expansion around the 
singularity and is done term-by-term. 

4. Weight-homogeneous vector fields 

We begin by introducing some useful notation and terminology. We write 

x(r) = ar p , (3) 

to denote the function 

x(t) = (ot* /St^V"), 

where a = (a, (3, 7, • • •) € C n \ {0}, and p = (p, q, r, ■ ■ •) G Q ra . A function of the form © is scale 
invariant in the sense that a change in the time scale, r' = nr, reveals that we must also 

have x' = r/ p x, and conversely. 

Now we demand that a scale invariant function x(r; xo) of the form (j3J) is an integral curve of 
the vector field f passing through xo , or equivalently, the associated dynamical system x = f (x) 
has a particular solution that is scale invariant, valid in a neighborhood of the assumed finite- 
time singularity. This means that 

apr p 1 = f(ar p ). (4) 

The notation ap stands for the monomial aiPi and is valid for each i = 1, • • • , n (range convention, 
no summation). 

We say that a vector field f (or the associated dynamical system) is scale invariant if it 
satisfies 

f(ar p ) = T p_1 f(a). (5) 

More generally, a vector field is called weight-homogeneous with weighted degree d if there is a 
vector p £ C n \ {0}, called the weight, and a vector d = (d\, ■ ■ ■ , d n ) such that 

f( ar p ) = r d f(a). (6) 

We denote the degree by deg(f , p) = d. When p = 1 the vector field is called homogeneous of 
degree d. Note that any scale invariant vector field is a weight-homogeneous field f = (fx, ■ ■ ■ , f n ) 
such that each component fi has degree pi — 1, from which it follows that dfi/dxj has weight 
Pi — 1 — pj. If f has degree d, then df/dxj has degree d — Pjej, with e.,- the j-th unit vector. 
There are weight-homogeneous vector fields that are not scale invariant; for instance, all linear 
dynamical systems with constant coefficients which are homogeneous vector fields. 

Using Eqs. (@J), © we see that given any nonzero vector p, a scale-invariant vector field f 
admits a scale-invariant integral curve provided that the inhomogeneous linear system 

ap = f(a), (7) 

has nontrivial solutions for a. Note that there may be nontrivial solutions for a with some 
components zero. When at least one of the components of a is nonzero, the corresponding 
solution of the form (j3J) is a particular solution of the dynamical system. Therefore, in this 
case, we know the exact asymptotic behaviour of the vector field f in the neighborhood of a 
finite-time singularity, given by the solution (j3J) with suitable a's and p's. 



5. Vector-field decompositions 

Unfortunately, most vector fields are neither scale invariant nor weight-homogeneous. However, 
since any analytic vector field f(x) can be expanded in a power series in some domain T> C M n , 
by taking any xo £ T> for any x € T> distinct from xo, any such vector field can be decomposed 
into weight-homogeneous components by taking for instance the first k + 1 terms in its Taylor 
expansion around xo- A simpler example of a vector field admitting a weight-homogeneous 
splitting is to take f to be any polynomial vector field; then, a possible decomposition is to split 
it such that each component is a suitable combination of monomials. 

We say that the nonlinear vector field f on Ai n admits a weight-homogeneous decomposition 
with respect to a given vector p if it splits as a combination of the form 

f = f (o) +f (i) + ... + f « (8) 

where the components f v\ j = 0, ■ • • , k, are weight-homogeneous vector fields, namely, 

fO')(ar p ) = T P+ 1 (9°' ) - 1 )fO')(a), j = 0, • ■ • , k, (9) 

for some non-negative numbers q^) and all a in some domain £ of W 1 . In terms of individual 
components, condition Q reads 

^ ) (arP) = T^ (j) - 1 /i J) (a), i = l,---,n, j = 0,.--,k. (10) 

In a slightly vague but suggestive manner we can say that a weight-homogeneous decomposition 
splits the original vector field in parts starting by collecting together the most nonlinear part 
and then proceeding down to the 'weakest' component such that each term in the splitting is 
'less nonlinear' than the previous one in a precise sense. 

There are two important features of such a vector- field decomposition (JHJ), @. Firstly, it 
is not unique. In general, for the given vector field f , many different vectors p can be found 
which each lead to different weight-homogeneous decompositions of the vector field in the same 
domain V. Secondly, since the subdominant exponents j = 1, ■ ■ ■ , k can be ordered, 

= g(°) < g &'0 < q^ 2 \ when j x < j 2 , (11) 

the degrees of the component vector fields in the decomposition Q are also ordered and (only) 
the first vector field f (°) appearing in the decomposition is scale invariant. 

Therefore, a weight-homogeneous decomposition of a vector field with respect to a vector p is a 
splitting into k+l weight-homogeneous components each with degree p+l(gk>) — 1), j = 0, • • ■ , k, 
such that the lowest-order vector field in the decomposition is scale invariant. The lowest-order 
vector field, f in a splitting is sometimes called the dominant part of f and includes the 
most nonlinear terms in f , whereas the remaining sum of parts, f sub = £j=i f {j \ is called the 
subdominant part. 

Given a vector field f , it is very important to have the complete list of all possible weight- 
homogeneous decompositions it admits; in other words, to know all the possible dominant and 
subdominant ways it can be split. The asymptotic method we employ to trace the behaviour 
of vector fields and their integral curves in a neighborhood of a movable finite-time singularity 
(of a particular or general solution), begins by finding all weight-homogeneous decompositions 
of the vector field valid in that neighborhood. 

6. The dominant balance 

Suppose that there exists a decomposition 



f = f(o) +f (D + ... + f « 



(12) 



into (k+1) weight-homogeneous components, such that the dominant part f is scale invariant, 
and each subdominant component f ^0, j = 1, • • • , k is weight-homogeneous. The scale invariant 
solution 

x (0) (r)=ar p , a ^ 0, p G Q", (13) 

is a solution 3 of the dominant part x = f(°) of the vector field provided that ap = f"(a). 
Thus, some of the components of the vector a may vanish. We call the components of the 
vector p = (pi, ■ ■ ■ ,p n ) the dominant exponents. In this case, we sometimes say that (|13j) is an 
asymptotic solution of the original system x = f(x). 

On the other hand, the k subdominant components f ^) satisfy 

/ i (i) (arP)=r^ U) - 1 /?' ) (a), i = l,---,n, j = l,--,k, (14) 

with the subdominant exponents j = 1, which are ordered and strictly positive. 

Dividing both sides by r p_1 and taking the limit as r — > 0, the true meaning of the subdominant 
exponents is revealed, and we have 

l im t ^Li = o, (15) 

which proves that the subdominant part of the vector field, f su , is less dominant than the 
dominant part, f^°\ which of course asymptotes as r p_1 . 

We say that the pair (a, p) is a (dominant) balance for the vector field f if the latter admits 
a decomposition satisfying Eqs. Q12JI - Q14JI . There may be several different balances for any 
particular decomposition of f and a way to classify them is by means of their order. The order 
of a balance (a, p) is the number of the nonzero components of the vector a. For a vector field on 
A4 n , the highest order of any possible balance is n and in this case the scale-invariant solution 
(|13|) corresponds to a possible dominant behaviour of a general solution of the original system 
x = f near the singularity. On the other hand, balances of a lower order than n describe possible 
asymptotics of particular solutions. 

There is an elegant convex-geometric explanation of the dominant balances of a vector field 
f . This requires us first to express the vector field in the so-called quasi-monomial form. Then, 
a dominant balance of order d corresponds precisely to a d-dimensional face of the Newton- 
Puiseux-Bruno polyhedron associated to f (cf. [0], [Zj). 



7. The K-matrix 

We can now study an important square matrix, the so-called (K-)ovalevskaya matrix, associated 
with a given vector field f. Consider the dominant part, f^, of f which admits an exact solution 
of the form (|13|) . described by the dominant balance (a, p) (of any nonzero order). The K-matrix 
of the vector field f at the balance (a, p) is the square matrix 

/C = Df (0) (a)-diagp. (16) 

The (K-)ovalevskaya exponents associated with the balance (a, p) are the n eigenvalues 
' • • )Pn) of /C. When the order of the balance is n (aj ^ for all i = 1, • • • , n), the K- 
exponents are called the resonances of /C. 

Setting w = f(°)(ar p ), and differentiating with respect to r, we have 

W = (apT^ 1 )- = diag(p - l)apr p ~ 2 , (17) 

3 due to a theorem of Goriely and Hyde U, 0, we can restrict attention to real vectors a only, as complex ones 
do not describe the behaviour near finite time singularities in the sense of Eq. 



while by the chain rule 



w = Df^arPjaprP- 1 

= DfW(a)a P T p - 2 , (18) 

where the last equality is most easily understood if we expand the derivative to obtain 
Wi = r Pi ~ 2 V//°^(a) ■ ap, for the i-th component. Thus, from the last two equations, we find 

/Cap = — ap; (19) 

that is, the K-matrix always has the vector f(°)(a) as an eigenvector with eigenvalue equal to 
pi = —1. We say that a balance is hyperbolic if the remaining (n — 1) K-exponents have positive 
real parts. 

Suppose now that we know all the eigenvectors and eigenvalues (— l,p2, • • • ,pn) of the 
K-matrix. By simply inspecting the form of the two sides in the variational equation for the 
dominant part of the vector field, namely, the equation 

w = L>f (0) (aT p )w, (20) 

we can write the set of fundamental solutions w'*' of this linear equation for w in the form 

W W = g W (log t)t p+ P\ (21) 

where, depending on whether or not K, is semi-simple, the gW(logr)'s are the eigenvectors 

or, in general, polynomials in logr. Hence, the solutions of the variational equation will be 

appropriate sums of terms by the form 1)21(1 . 

Therefore, any solution of the original system will be well approximated (cf. 8 , p. 299, for 
the precise conditions) by a solution of the form (this is a Taylor estimate): 

x = x (0) + w, (22) 

where 

k 

x (°) (r) = ar p , w = ^ hW r P+« , (23) 
i=i 

where the h^'s are, in general, polynomials in logr. Furthermore, we arrive at the interesting 
conclusion that in any particular or general solution, the arbitrary constants characterizing it 
will first appear in those terms whose coefficients have indices equal to a K-exponent. In a 
general solution, the arbitrary constants normally appear at different places in an expansion, 
and consequently, a solution in which a K-exponent (different from the —1 value which as we 
have shown always exists) is either negative, or has nontrivial multiplicity, may or may not 
be a general solution. We shall see in the next section how we can use Eq. (|22jl to obtain 
an important series representation of the solutions to the original dynamical system near a 
finite-time singularity. 

Apart from the principal use made here in unravelling the nature of finite-time singularities, 
there are various important connections between the K-exponents, first integrals, integrability 
properties of Hamiltonian systems and complex algebraic geometry, cf. |10| . 



8. Formal expansions 

In the previous subsection, we derived the formal expansion (1221) for the solution of the dynamical 
system x = f (x) near a finite-time singularity by including only the first terms in a power-series 
expansion of x. That solution corresponds to a given dominant balance T = (a, p) and depends 
on the K-exponents (— 1, p%, ■ ■ ■ , p n ) associated with the dominant part, f(°\ of the original 
vector field f. Writing Eq. 1)22(1 in full, that is including terms of all orders, can lead to very 
complicated expansions in general, and the essential feature is the appearance of logarithmic 
terms. We note that there are a number of general theorems guaranteeing the existence and 
convergence of such expansions - see, for instance, jllj for a review. 

Let £ + be the set of all K-exponents with positive real parts. For simplicity we assume that 
all K-exponents in £ + are rational, and define the number 1/s to be the least common multiple 
of the denominators of the numbers in the set Ti = {q^\ ■ ■ ■ , q^} U £ + . In the case where 
any log terms are absent (for instance, when K is semi-simple), we can write, by (|22|). the full 
expansion of the general solution around the finite-time singularity in the form of a Puiseux 
series, 

x = r p U + £c i r i /M, (24) 

where, as we know already from the previous section, each of the n arbitrary constants in (|24[) 
will first appear in the term with coefficient c^, k = ps and p € £ + . 

Hence, finding the final form of the solution (general or particular, depending on the number 
of arbitrary constants appearing in the series expansion) can be now reduced to knowing the 
coefficients Cj in the expansion. These coefficients are computed by inserting the Puiseux series 
(|24[) into the original system. This leads to a set of recursion relations, a linear system for the 
coefficients Cj. For the jth-order coefficient we find 

j 

JCcj - -Cj = Pj(ci,- • ',Cj-x), (25) 

where the forms Pj are polynomial in its variable, read off from the original equation. 

There is an important consistency condition to be satisfied for the above analysis to be valid. 
Multiplying both sides of ()25|) by v, an eigenvector of the K-matrix, we see that when j/s = p, 
an eigenvalue of the K-matrix, we must have the following compatibility condition (v T denotes 
the adjoint eigenvector of v): 

v T -P j = 0, forallpe£ + . (26) 

Therefore, if the above compatibility condition is violated at some eigenvalue, then we conclude 
that no solution in the form of a Puiseux series can exist and we have to search for more general 
solutions which may contain logarithmic terms. Such a more general series will be of the form 
of a ^-series (cf. j!2j for this terminology): a direct generalization of the form (|24[) 



x 



+ EE c ^ /s ( rP1 °g r ) J/s l' (27) 

i=l 3=1 



where p is the first K-exponent for which the compatibility condition is not satisfied and s as 
defined above. The procedure for the calculation of the coefficients in this more general case is 
the same as before, leading again to the form of the general solution in a suitable neighborhood 
of the finite-time singularity as a ■(/'-series. 



9. Relation to unstable manifolds 

Suppose that (|T3j) is an asymptotic solution of the vector field x = f (x). We set 



x(t)=tPX(s), S = logr, (28) 

and imagine that an equation of the form 

X(s) = const. = a, (29) 

regards the coefficient a of a given balance (a, p) of the vector field f as an equilibrium point of 
the new system given, in terms of the new variables X(s) = (Xi(s), • • • ,X n (s)), by 

X[ = F l (X u ---,X n ), i = l,...,n 
X' n+1 = qX n+1 , X n+1 =ei s , (30) 

where (') = d/ds and q is the least common multiple of the denominators of the subdominant 
exponents in the original system. We call the dynamical system (|3()[1 the companion system 
of the original vector field x = f(x). Thus, the transformation (|3()j) to the companion system 
associates a different system of this sort to each one of the balances (a, p) of the original system 
provided that any given pair of balances has different dominant exponents p). 
Consider now the linearized system 

X' = DF(a)X, (31) 

which is the variational equation of the companion system around the equilibrium point 
X = (a, 0). This is a constant-coefficient, linear system. From the fundamental theorem of such 
systems, it follows that the general solution passing through the initial condition X(0) = Xo is 
given by 

X(s) = e sDF(a) X . (32) 
Now we know (cf. ^lj) that for any n x n matrix A: 

e sA v = J2^ 8 Vj(s), (33) 

where Xj are the eigenvalues of A and Vj are C n -valued polynomials, the latter being constant 
if and only if A is semi-simple (diagonalizable) . So 

X = J>,-(log t)t x i, (34) 

in terms of the r-time. Hence, using this form of the local solution to the companion system, 
we can express the local solution to the original system around its finite-time singularity in the 
form 

x W = E^ lo s r ) rP+Aj > ( 35 ) 

which is precisely the Taylor estimate leading eventually to the psi-series representation we 
arrived at in previous Sections. We therefore reach the interesting conclusion that around the 
equilibrium point (0, q), origin of the companion system the eigenvalues of that system are 
simply (p, q); that is, the dominant exponents of the asymptotic solution of the original system 
together with the number q characterizing the subdominant part of the vector field. Moreover, 
around the equilibrium point (a, q) of the companion system (|3T)|) . the eigenvalues Ai, • • • , Xj,q 
of the companion system are (pi, ■ ■ ■ ,p n ,q); that is, they are precisely the K-exponents of the 
original system together with the subdominant number q. 



We have an interpretation of the results of the previous Sections concerning the local 
behaviour of the original system around its finite-time singularities from a dynamical systems 
perspective. Using the companion transformation, a local analysis of the companion system 
X' = F(X) around its equilibrium points (with s — » — oo necessarily) will provide the local 
analysis of the solutions of all possible balances of the original system around its singularities 
(r = e qs — > since q > always, for an acceptable decomposition). 

Note that, since we need s — > — oo, we are only interested in the unstable manifold (eigenvalues 
with positive real parts) of the equilibrium points of the companion system. Therefore, the 
negative K-exponents (corresponding to the stable manifold of the companion system) are not 
connected to the behaviour of solutions of the original system at the finite-time singularity, but 
are associated with its behaviour as r — > oo. 

For a discussion of the companion transformation in connection with integrability and 
complex dynamics, see and Refs. therein. 

10. The method of asymptotic splittings 

The results in the previous subsections suggest a general procedure to uncover the nature of 
singularities by constructing series expansion representations of particular or general solutions 
of dynamical systems in suitable neighborhoods of their finite-time singularities. This method 
consists of building splittings of vector fields that are valid asymptotically and trace the dominant 
behaviour of the vector field near the singularity. A resulting series expansion connected to a 
particular dominant balance helps to decide whether or not the arrived solution is a general one 
and to spot the exact positions of the arbitrary constants as well as their role in deciding about 
the nature of the time singularity. 

The method we suggest below is analogous to the so-called ARS procedure connected with 
the Painleve and integrability properties of dynamical systems |13j . but here the whole approach 
and viewpoint are completely different. We are not concerned with notions of integrability but 
solely with the problem of the nature of finite-time singularities. Also, our systems are real- 
valued with a real time variable. 

To apply this method of asymptotic splittings to a particular dynamical system in an effort 
to discover the nature of its time singularities, we must follow this recipe: 

(i) Write the system of equations in the form of a dynamical system x = f(x) with x = 
(xi, • • • , x n ), and identify the vector field f (x) = (/i(x), • • • , / n (x)). 

(ii) Find all the different weight-homogeneous decompositions of the system; that is, the 
splittings of the form 



in 



) 






into the equation x = f(°). Study the resulting algebraic systems, and find all dominant 
balances (a, p) together with their orders. 

Identify the non-dominant exponents, that is the positive numbers qv\ j = 1, • • • , k, such 
that 




f (0) _> £>f (0) 



£>f<°>(a) -» £>f (0) (a) - diagp. 



(vi) Compute the spectrum of /C, 



spec(7C) = (-1,P2, • • • , Pul- 
ls K, semi-simple? Are the balances hyperbolic? 

(vii) Find the eigenvectors of IC. 

(viii) Identify s as the multiplicative inverse of the least common multiple of all the subdominant 
exponents and positive K-exponents. 

(ix) Substitute the Puiseux series 

oo 

Xi = ^2 Cji T Pt+ i 
3=0 

into the original system. 

(x) Identify the polynomials Pj and solve for the final recursion relations which give the 
unknown coefficients Cj. 

(xi) Check the compatibility conditions at the K-exponents, 

vj ■ P p = 0, for each eigenvalue p. 

(xii) If the Puiseux series is valid, then the method is concluded for this particular splitting. 
Otherwise, if compatibility conditions are violated at the eigenvalue p*, restart from step 
ED by substituting the logarithmic series (|2TJ) . 

(xiii) Get coefficient at order p*. Write down the final expansion with terms up to order p* . 

(xiv) Verify that compatibility at p* is now satisfied. 

(xv) Repeat whole procedure for each of the other possible decompositions. 

Note that although the whole spectrum of possible behaviours of the system near a time 
singularity is concluded once we find valid series expansions corresponding to each balance 
in each particular decomposition, an additional analysis of the phase space of the companion 
systems corresponding to each one of the balances of the original system may lead to valuable 
insights as to the geometric structure of the phase space - how the orbits behave - in the 
neighborhood of the finite-time singularity. 

Following the above steps even up to that of calculating a dominant balance in one particular 
decomposition, can be very useful since it gives you one particular possible asymptotic behaviour 
of the system near the time singularity. In this respect, the whole method expounded here is 
truly generic since it helps to decide the generality of any behaviour found in an exact solution - 
that is, how many arbitrary constants there are in the final solution that shares that behaviour 
(particular or general solution). It is rare that a Puiseux series is inadequate to describe the 
dynamics (semi-simplicity of K), but in such uncommon cases one must resort to the more 
complex logarithmic solutions. 



11. A worked example 

As a concrete example of the above analysis, consider the homogeneous and isotropic FRW 
cosmological equations in general relativity for a perfect fluid of pressure p(t), density p(t), and 
equation of state p = wp, where w is a constant. They read 

-~(l + 3w), (36) 
0, (37) 
p-4 (38) 



a 
a 

p + 3Hp(l + w) = 
3H 2 = 



where a(t) is the scale factor, H(t) the Hubble expansion rate and k the constant spatial 
curvature. Setting a = x,x = y, p = z, this system reads 



x = y, (39) 

y = axz, (40) 

z = byz/x, (41) 

with 

a = -(1 + 3m)/6, b = —3(1 + «;), (42) 
and is subject to the integral constraint (the Friedmann equation) 

3^ + \ = z. (43) 

This system is weight-homogeneous. The unique balance is determined by 

/ 2(2 + b) \ (11 \ 

a = (a, ft 7) = ( a, pa, ^ J , p = (p,g,r) =(--,--- 1,-2 J , (44) 

and the parameters a, 6 depend on via Eq. (]42[). Since the field is weight-homogeneous (i.e., 
f = f(°)), this balance corresponds to an exact, scale-invariant solution of the original system. 

Note that there is one arbitrary coefficient, and so we expect that one of the K-exponents will 
be zero. For the vector field f = (y, axz, byz/x) the associated K-matrix, K = Df(a) — diagp, 
has characteristic equation 

r 3 + (2p - l)r 2 + 2(p - l)r = 0, (45) 

hence the K-exponents are 

r = -1, 0, 2(1 -p). (46) 

We can use the integral constraint 1)43(1 to get a value for the coefficient a. Balancing the terms 
in (@3J) leads to 

In the case when all components of the vector a are real, the solution X = (x, y, z) = (a, a, p) of 
the dynamical system experiences a finite-time singularity. When 1 — p > 0, that is when either 
w < —1 orw > —1/3, solutions are general; while, in the range —1 < w < —1/3, we have only 
behaviours corresponding to particular solutions of the system. Note that when w = 0, we have 
b = —3 and we find a behaviour similar to that of dust, but when w = 1/3, b = —4 the behaviour 
is that of the standard radiation models. Further, calculating the recursion relations to compute 
the coefficients of the series expansion term by term, we find the following asymptotic solution 
for radiation, 

X = arV 2 _ | T 3/2 _ ^ C3lT 5/2 + . . . 

while the dust-dominated expansion is found to be 

x = ar 2/3_« r 4/3 + ... 

4 



4 



2 i ^-4/3 



The method that we have described in detail here provides a toolkit for the investigation 
of the general form of a range of finite-time singularities in general relativistic cosmologies. In 
particular, the 'sudden' singularities introduced by one of us JH]> in which a, a, and p remain 
finite but a and p — > oo at finite time in situations where no functional relationship is assumed 
between p and p, have been widely studied ^2]- Such situations, allow sudden singularities to 
develop at finite time without violating the strong-energy condition of general relativity and are 
require less severe conditions than future 'big rip' singularities |17j with w < — 1. Elsewhere, 
we will report on the results of applying these methods to determine the general behaviour on 
approach to these singularities in general relativity and in higher-order gravity theories. 
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